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SUMMARY 


Smoke  visualization  and  embedded  laser  velocimeter  operating  in  the  backscatter  mode  were  used  to 
survey  the  flow  about  wings  oscillating  through  stall  in  the  Luminy  wind  tunnels  S1L  and  S2L. 

The  experiments  were  carried  out  about  rectangular  wings  with  an  NACA  0012  airfoil  section  at 
nominal  Reynolds  numbers  R^IO5,  2105,  and  106.  During  the  oscillation,  the  wing  angle  of  attack  varied  around 
mean  values  oto=6deg  and  15deg  through  an  angle  amplitude  Aa=±  6deg,  at  reduced  frequency  k=0.188. 

Pictures  and  movies  realized  using  the  laser  sheet  technique  have  shown  that  the  features  of  dynamic 
stall  are  independent  of  Reynolds  numbers,  at  least  in  the  present  experiments  (105<R<j<106).  Effects  of  upstroke 
and  downstroke  on  separation  and  reattachment  of  boundary  layer  are  clearly  shown  as  well  as  deviation  from 
steady  behavior  analyzed  at  the  same  angle  of  attack  for  fixed  and  oscillating  airfoil. 

The  measurements  of  two  orthogonal  components  of  instantaneous  velocity  through  the  boundary 
layer  along  a  normal  to  the  surface  of  the  wing  have  shown  that  the  embedded  laser  velocimeter  is  well  suited  to 
characterize  qualitatively  and  quantitatively  the  structure  of  the  boundary  layer  (laminar,  transition,  turbulent),  as 
well  as  strong  separation  with  leading  edge  bubble  involved  during  dynamic  stall.  These  results  constitute  a 
useful  data  base  for  future  computational  unsteady  boundary  layer. 

List  of  Keywords:  Dynamic  stall,  Boundary  layer,  Bubble  separation,  Laser  light  sheet  visualization, 
Embedded  Laser  Doppler  Velocimeter,  Image  processing. 


INTRODUCTION 

The  loading  of  an  airfoil  during  dynamic  stall  has  been  extensively  investigated.  The  features  of  the 
stalling  process,  typical  of  the  retreating  blade  of  a  helicopter  rotor  in  forward  flight,  have  been  delineated 
several  years  ago  on  airfoils  oscillating  either  in  pitch(1'3)  or  in  fore  and  aft-pitching  combined  motion(4'6).  More 
recently,  W.  G.  Bouseman(7)  has  shown  that  the  loading,  characterized  by  the  peak  airfoil  lift  and  drag,  and  the 
minimum  pitching  moment,  termed  the  “dynamic  stall  function”,  is  similar  over  a  wide  range  of  test  conditions 
including  flight  tests.  In  particular,  the  dynamic  stall  function  appears  to  be  relatively  insensitive  to  Reynolds 
number  and  shows  that  an  accurate  prediction  of  3-D  effects  is  not  of  prime  importance  for  the  prediction  of 

dynamic  stall  on  a  flight  vehicle. 

These  results  encourage  to  enforce  numerical  methods  presently  developed  for  the  direct  calculation 
of  dynamic  stall  by  solving  the  Navier-Stokes  equation(9)  but  limited  to  2D  flow  at  low  Reynolds  Number.  These 
challenging  methods  are  at  a  research  stage  and  their  validation  and  improvement  require  more  accurate  and 
suited  experimental  data.  Particularly,  phenomena  occurring  in  the  region  very  close  to  the  blades,  like  bubble 
formation,  transition  and  separation,  require  local  and  detailed  information  from  the  boundary  layer  (BL). 
Among  these  informations,  visualization  and  instantaneous  velocity  profiles  remain  one  of  the  main  useful  and 
suitable  tools  for  analyzing  the  behavior  and  the  structure  of  BL. 

The  present  experimental  study,  done  at  LABM  (Laboratoire  d’Aerodynamique  et  de  Biom^canique 
du  Mouvement)  from  the  University  de  ia  Mediterranee,  proposes  to  contribute  to  the  numerical  and  theoretical 
developments  undertaken  in  this  research  field  at  the  Aeromechanics  of  the  Army/NASA  Rotorcraft  Division. 
The  experimental  results  provided  on  BL  velocity  profiles  and  on  the  structure  of  unsteady  separated  flows  have 
been  mainly  based  on  measurements  using  a  unique  embedded  LDV  technique.  LABM  has  developed  for  5 
years  original  types  of  embedded  LDV  for  measurements  of  oscillating  profile  and  rotary  wings  boundary 
layer(6),(9\  Concerning  oscillating  profiles(6)  (NACA0012,  0.3m  in  chord),  the  velocity  components  have  been 
measured  in  2D  and  3D  configurations  by  help  of  a  LDV  embedded  in  the  oscillating  frame,  in  such  a  way  that 
the  velocities  are  directly  measured  in  a  reference  frame  in  relative  motion.  Concerning  hovering  helicopter  rotor 
blades(9\  BL  measurements  have  been  performed  in  a  frame  linked  to  the  rotating  blade  for  one  radial  distance 
from  the  rotation  axis  and  at  a  chord  abscissa  x/c=0.25. 

The  aim  of  the  present  seed  proposal  (first  phase)  was  to  show  that  the  embedded  LDV  applied  to  a 
blade  profile  oscillating  through  dynamic  stall  is  well  suited  to  study  qualitatively  and  quantitatively  the  physics 
of  unsteady  separated  flows. 

This  preliminary  study  has  been  undertaken  in  2  different  wind  tunnels  (S1L  and  S2L  of  LABM)  at  3 
Reynolds  numbers  (105,  2.105  and  106)  corresponding  to  3  Mach  numbers  (0.015,  0.030  and  0.15).  The  airfoil 
was  oscillating  in  pitch  through  dynamic  stall  around  a  mean  angle  of  attack  Oo,  at  a  reduced  frequency  k  and  an 
amplitude  of  incidence  Aa  relevant  to  the  retreating  helicopter  rotor  blade  (k=0.188,  ao=6deg  and  15  deg, 
Aa~6deg). 

A  particular  attention  has  been  paid  in  this  first  phase  to  the  structure  of  the  flow  field  when  the  angle 
of  attack  is  increasing  and  decreasing  through  the  separated  and  reattached  BL  regimes. 


The  formation  and  shedding  of  vortices  generated  by  the  dynamic  stall  have  been  observed  using  a 
laser  sheet  visualization  technique.  Movies  showing  separation  and  reattachment  during  dynamic  stall  have  been 
compared  at  Re=105  2.10s  and  106. 

The  two  velocity  components  (U=U(<Dt),  and  V=V(fi)t))  have  been  measured  in  the  BL  from  y=0.3mm 
at  the  wall  to  100  mm  outside  of  the  BL.  Inside  the  BL,  measurements  have  been  performed  at  every  altitude 
step  of  0.2  mm. 

SYMBOLS 

C  Airfoil  chord,  m 

f  frequency,  Hz 

h  Airfoil  span,  m  > 

k  reduced  frequency,  k=Cco/2Uoo 

Re  Reynolds  number,  Re=UooC/v 

Uoo  Wind  tunnel  velocity,  m/sec 

U,V  Tangential  and  normal  velocity,  m/sec 

s  Curvilinear  abscissa,  m 

t  Time,  sec 

x,  y  Tangential  and  normal  coordinate  at  the  airfoil  surface,  m 

a  Angle  of  attack  of  the  airfoil,  deg 

a0  Mean  angle  of  attack,  deg 

Aa  Amplitude  of  the  pitching  oscillation,  deg 

8  Thickness  boundary  layer,  m 

r|  Reduced  normal  coordinate,  r|=y(Rex)1/2/x 

©  Rotational  frequency,  rad/sec 


EXPERIMENTAL  INVESTIGATION 
The  wind  tunnels 

-SILuminy  high  speed  (Fig.l)  is  a  closed  circuit  wind-tunnel  with  octagonal  closed  test-section 
(inside  circle  of  3m),  6m  in  length,  free  stream  velocity  varying  from  5m/s  to  lOOm/s;  constant  temperature 
of  the  flow  within  the  test  section  :  At<l°C;  natural  turbulence  intensity  lower  than  0.3%  (longitudinal)  and 
0.1%  (transversal). 

-S2Luminy  low  speed  (Fig.2)  is  an  open  circuit  wind-tunnel  with  rectangular  closed  section 
(lmxO.Sm)  3m  in  length,  free  stream  velocity  varying  from  2.5m/s  to  25m/s,  natural  turbulence  intensity  less 
than  0.5%. 


Oscillating  airfoils 

Airfoils  used  in  both  wind  tunnels  are  of  NACA0012  profile.  One  of  them,  experimented  in  S1L,  is 
0.488m  in  chord,  1  m  in  height,  with  one  tip  mounted  flush  to  the  floor  of  the  test  section,  the  other  tip  is 
equipped  with  an  end  plate  in  order  to  simulate  a  like  2-D  flow  at  mid  span.  The  pitching  motion  of  the 
airfoil  is  allowed  by  means  of  an  oscillating  device  located  under  the  test  section  (see  Fig.l)  and  supporting 
vertically  the  model  by  means  of  a  support  shaft  located  at  the  quarter  chord  axis.  In  this  way,  the 
instantaneous  incidence  a(t)  of  the  airfoil  is  given  by: 

a(t)=  a0+Aasin((Dt), 

where  Oo  is  the  mean  incidence  and  g>  the  rotational  frequency. 

The  airfoil  used  in  S2L  (0.30m  in  chord  and  0.50m  in  height)  is  mounted  vertically  and  spans  the 
entire  test  section.  This  airfoil  is  supported  at  the  quarter  chord  axis  by  an  oscillating  device  (see  Fig.2) 
which  allows  a  combined  motion  of  fore  and  aft  and  pitching  motion.  In  the  present  case  only  the  pitching 
motion  was  used,  giving  the  same  variation  of  the  instantaneous  incidence  as  in  S1L. 

The  following  table  summarizes  the  experimental  conditions  realized  in  S1L  and  S2L: 


Table  1 


C 

8 

1 

Chord  C:  m 

Re 

f=27i/o:  Hz 

k 

a0:  deg 

Aa:  deg 

S1L 

31 

0.48 

io5 

3.8 

0.188 

21 

6 

S2L 

5 

0.3 

lo’^io" 

1 

0.188 

6,  15 

6 

Visualization 

The  technique  consists  in  illuminating  the  mid:span  of  the  airfoil  while  the  flow  is  seeded  upstream 
with  oil  smoke  emitted  through  a  tube  rake.  The  light  sheet  is  produced  from  an  optical  system  mounted  flush 


with  the  wall  of  the  wind  tunnel  in  order  to  avoid  light  reflection.  The  optical  system,  manufactured  by 
Dantec,  is  composed  of  a  glass  cylinder  and  a  converging  lens  connected  by  an  optical  fiber  to  a  laser  source 
(Spectra-Physics  12  Watts).  Figure  3  shows,  as  an  example,  the  set-up  realized  in  S2L.  It  can  be  seen  that  the 
light  sheet  has  been  adjusted  to  illuminate  the  middle  horizontal  plane  of  the  wind  tunnel.  During  tests, 
pictures  of  the  flowfield  close  to  the  airfoil  are  taken  from  a  digital  camera  (Olympus  Camedia  C-2500L) 
synchronized  on  the  motion  of  the  oscillating  airfoil.  These  pictures  are  processed  through  “Adobe 
Photoshop44  software  in  order  to  capture  the  features  of  the  flowfield.  Moreover,  movies  have  been  realized 
from  the  pictures  shot  at  different  instant  of  the  period  of  oscillation. 

The  same  technique  described  above  has  been  used  in  S1L. 

Velocity  measurements 

For  each  chordwise  section  selected,  the  flow  survey  is  performed  along  the  local  normal  to  the  wall 
surface  (the  altitude  Y  ranging  from  about  0.3  mm  to  100mm). 

The  Embedded  Laser  Doppler  Velocimeter(6)  (ELDV)  used  for  this  survey  has  an  optical  head 
mounted  on  a  supporting  turntable  linked  to  the  oscillating  frame  as  sketched  in  Fig.l  and  2.  Moreover,  it  is 
equipped  with  a  beam-expander  to  increase  the  focal  distance  to  400mm.  This  optical  head  is  installed  on  an 
automated  2D-displacement  device  mounted  itself  on  the  turntable.  The  laser  beams  are  so  focusing  from  a 
45  deg  mirror  in  the  boundary  layer  at  the  span  wise  locations  considered.  The  supporting  turntable  is  linked 
with  the  oscillating  frame,  so  that  U  and  V  velocity  components  can  be  directly  measured  in  the  same 
reference  frame  than  the  oscillating  BL.  A  teledriven  system  allows  the  adequate  positionning  of  the 
measurement  volume  at  any  point  of  the  airfoil  surface  (30cm  in  chordwise  displacement).  An  angular  sector 
provides  the  selection  of  the  surveying  normal  direction,  and  the  laser  measurement  volume  can  be  displaced 
along  the  local  normal  to  the  surface  from  y  =  0.3mm  to  y  =  145mm  with  a  displacement  accuracy  of  0. 1mm. 
Due  to  the  periodicity  of  the  flow,  each  period  is  considered  as  a  specific  sample  of  the  same  phenomenon. 
So,  each  velocity  component  is  recorded  at  each  phase  angle  cot  comprised  between  Odeg  and  360deg  by  step 
of  ldeg  over  a  large  number  of  periods.  Data  are  then  statistically  analyzed  at  prescribed  values  of  the  period, 
e.g.  the  instantaneous  incidence,  with  a  precision  of  5a=4Aa/36(>= 24/360=0. 066deg.  The  acquisition  time  is 
stopped  when  19,564  data  have  been  stored,  that  generally  requires  about  150  periods. 

Data  acquisitions  are  made  on  a  microcomputer  from  two  Burst  Spectrum  Analyser  delivering  the 
values  of  the  two  components  U  and  V  and  the  arrival  validation  time  for  each  frequency  measurement.  The 
software  used  for  acquisition  and  data  reduction  (COMBSA)  has  been  developed  at  LABM  under  Apple 

LabVIEW  system. 

RESULTS  AND  DISCUSSION 
Visualizations 

The  visualizations  have  been  performed  owing  to  the  technique  described  above.  The  experiments 
have  been  duplicated  in  S1L  and  S2L,  in  order  to  investigate  Reynolds  numbers  ranging  from  105  to  106. 

Both  steady  and  unsteady  case  have  been  studied.  In  the  steady  case,  pictures  of  the  flow  on  the  upper 
side  of  the  profile  were  taken  from  the  digital  camera  (Olympus  Camedia  C-2500L)  for  values  of  the 
incidence  varying  from  Odeg  to  21  deg  by  step  of  ldeg.  Two  Reynolds  numbers  were  explored  (R*  =  10 5  and 
2.10  5)  in  S2L  and  one  (Rc  =106)  in  S1L.  Movies  representing  one  cycle  of  “quasi-steady”  oscillation,  has 
been  realized  from  the  successive  steady  pictures  corresponding  to  cto  increasing  from  a=  9deg  to  2 ldeg  and 
back  to  9deg  by  step  of  ldeg  in  S2L  and  from  15deg  to  27deg  and  back  in  S1L. 

In  the  unsteady  case,  pictures  of  the  upper  side  flow  field  were  shot  at  different  instant  of  the  period 
when  the  airfoil  oscillated  in  pitch  through  dynamic  stall  at  parameters  of  table  1,  and  movies  of  successive 
pictures  have  also  been  realized. 

Low  Reynolds  numbers 

As  examples,  Figures  4  to  7  shows  steady  and  unsteady  results  at  two  angles  of  attack.  One  of  them 
(a<f=10deg)  corresponds  to  unstalled  regime,  the  other  (a0=14deg)  to  stall  regime.  In  each  figure,  the  upper 
picture  corresponds  to  the  steady  case,  the  middle  one  to  the  pitching  airfoil  at  the  same  incidence  in  downstroke, 
the  lower  to  the  airfoil  in  upstroke.  Figures  4-5  are  relative  to  R^^IO5  and  Fig  6-7  to  R<.=2105.  Useful 
comparisons  can  be  made  from  these  figures:  for  the  range  of  variation  presented  here,  it  can  be  said  that  the 
effect  of  Reynolds  number  on  the  phenomenon  seems  not  to  be  significant.  Reattachment  of  boundary  layer 
during  downstroke  occurs  from  leading  edge  to  trailing  edge  and  is  very  close  to  steady  flow  (Figs  4  to  7), 
though  during  upstroke,  for  the  2  Reynolds  numbers  studied  (Fig.4  and  6),  separation  seems  to  occur  at  lOdeg 
with  the  formation  of  a  high  vorticity  flow,  sustaining  as  it  is  well  known  a  high  level  of  lift  until  the  maximum 
incidence  is  reached. 

Two  movies  showing  the  successive  phases  of  the  period  of  oscillation  have  been  realized:  one  when 
the  airfoil  oscillates  in  a  like  “quasi-steady”  flow,  the  other  when  the  airfoil  undergoes  the  unsteady  dynamic 


stall.  The  deviation  from  steady  behavior  due  to  unsteadiness  effect  involved  by  dynamic  stall  is  underlined  by 
comparison  between  the  two  movies. 

High  Reynolds  number 

Experiments  performed  in  S1L  have  concerned  a  rectangular  wing  of  the  same  NACA  0012  airfoile 
used  in  S2L.The  wing  of  0.86m  in  height,  and  0.48m  in  chord  spans  only  part  of  the  test  section.  It  seemed 

interesting  to  test  such  a  model  at  Re=10^  in  order  to  evaluate  both  the  effects  of  Reynolds  Number  and  3- 
dimensionnality  induced  by  the  wing  tip  by  comparing  the  results  to  those  obtained  in  S2L  in  2D  configuration 
at  lower  Reynolds  Numbers.  In  these  new  experimental  conditions,  static  stall  incidence  of  the  wing  is  about 
21  deg,  compared  to  12deg  in  2D  case.  The  parametric  conditions  of  pitching  around  a  mean  incidence  do 
-21 deg  are:  Aa=±6  and  k=0.188. 

Figures  8  toll  show  some  of  the  visualization  pictures  realized  at  incidences  varying  between  15deg 
and  27deg.  In  Figures  8  and  9,  the  effects  of  dynamic  stall  occurring  at  18deg  and  21  deg  during  downstroke  and 
upstroke  are  presented  and  compared  to  the  static  values.  It  can  be  clearly  seen  that  during  downstroke  the 
reattachment  of  the  boundary  layer  concerns  a  large  part  of  the  upperside  from  leading  edge  to  trailing  edge. 
During  upstroke,  a  separation  of  the  flow  consecutive  to  the  reattachment  process  is  characterized  by  a  high 
vorticity  level  generated  by  a  large  bubble  at  leading  edge  (a=21deg). 

In  Figure  10,  the  flow  is  reattached  at  lower  incidence  (a=16deg)  with  the  same  feature  than  in  steady 
configuration.  Nevertheless,  when  the  maximum  incidence  is  reached  at  27deg  (Fig.  11),  the  stall  experimented 
by  the  airfoil  in  unsteady  case  is  deeper  than  the  one  occurring  at  static  incidence  and  exhibits  the  high  vorticity 
flow  previously  mentionned. 

Visualizations  realized  at  different  Reynolds  Numbers  ranging  from  10^  to  10^  on  2D  and  3D  wings 
oscillating  at  the  same  reduced  frequency  reveal  the  same  features  of  the  dynamic  stall:  reattachment  during 
downstroke  from  leading  edge  to  trailing  edge,  and  separation  of  the  flow  with  strong  vorticity  during  the  end  of 
the  upstroke  generated  by  a  leadiung  edge  bubble. 

It  can  be  concluded  than  variation  of  Reynolds  effects  in  a  range  105-106  and  tip  wing  do  not  seem  to 
affect  the  structure  of  dynamic  stall  at  least  at  k=0. 188.  This  remark  is  in  good  agreement  with  some  conclusion 
given  in  the  paper  by  W.  Bouseman(7). 

Visualizations  have  shown  that  separated  flows  occurring  on  the  upper  side  of  the  airfoil  oscillation  at 
high  angle  of  attack  result  from  the  bursting  of  a  leading  edge  bubble.  On  the  other  hand,  when  the  airfoil 
oscillates  through  lower  incidences  with  no  separated  flow  (a<12deg),  a  bubble  remains  located  close  to  the 
leading  edge,  but  its  size  is  too  small  to  be  well  defined  by  the  image  processing.  However,  it  will  be  shown  here 
after  that  the  analysis  of  velocity  profiles  recorded  through  the  boundary  layer  can  result  in  the  quantitative 
definition  of  the  laminar  separation  bubble.  The  analysis  is  based  on  the  fact  that  when  the  flow  is  seeded 
upstream,  only  very  few  particles  can  enter  into  the  bubble.  Thus,  when  the  measuring  volume  is  inside  the 
bubble,  the  rate  of  data  acquisition  over  several  periods  remains  close  to  zero  resulting  in  a  no  data  zone 
observed  on  the  velocity  records.  This  no  data  zone  allows  determining  the  height  of  the  bubble  and  the 
incidence  at  which  the  bubble  enters  into  the  measuring  volume  location. 

Velocity  profiles  at  low  incidence  (0<a<12deg)  and  low  Reynolds  number  (Re=105) 

As  explained  before,  the  lack  of  data  on  time-velocity  records  reveals  the  presence  of  a  bubble. 
Velocity  components  U~U((ot)  and  V=V(cot),  have  been  recorded  for  different  location  Y  of  the  measuring 
volume  (Ymm-100;  80;  60;  40;  30;  20;  10;  8;  6;  4;  2;  1;  0.8;  0.6;  0.4),  at  different  chord  location  x/C=0. 1;  0.2; 
0.3;  and  0.67. 

Figure  12  gives  an  example  of  cumulative  (more  than  300  periods)  velocity  records  with  incidence  at 
x/C=0.1  obtained  at  different  distance  from  the  wall.  It  can  be  seen  no  lack  of  data  far  from  the  wall  (Y=10mm) 
when  the  incidence  varies  from  12deg  to  0  and  back  to  12.  When  approaching  the  wall  (from  Y=3mm  to 
Y=0.6mm),  a  zone  of  “no  data”  indicative  of  a  separation  bubble  can  be  observed  during  downstroke  and 
upstroke.  It  has  thus  been  possible  to  determine  the  height  of  the  bubble  at  a  given  incidence  for  different  chord 
locations  from  0. 1  to  0.67. 

Figure  13  shows  the  variation  with  the  airfoil  incidence  of  the  size  of  the  separation  bubble  during 
downstroke  and  upstroke  deduced  from  the  different  velocity  records.  It  can  be  observed  that  for  0deg<a<12deg, 
the  separation  bubble  remains  located  near  the  leading  edge  (0<x/C<0.3),  as  no  bubble  has  been  identified  at 
x/C=0.3  and  0.67  (see  Fig.  13).  The  maximum  size  of  the  bubble  can  be  evaluated  and  is  shown  to  vary  between 
6mm  (Y/C=0.02)  in  downstroke  and  3mm  (Y/00.01)  in  upstroke. 

Velocity  profiles  U(Y,a),  V(Y,a),  have  been  statistically  analyzed  from  the  time-velocity  records  at 
different  instants  of  the  period  corresponding  to  different  values  of  airfoil  incidence  a.  Although  useful 
information  could  have  been  carried  out  from  the  V(Y,a)  component  that  constitutes  with  U(Y,a)  a  valuable 
data  base  to  be  analyzed  in  more  detail  later,  only  some  typical  results  of  the  U(Y,a)  component  have  been 
analyzed  in  the  frame  of  this  report.  Figure  14  shows  as  an  example  the  results  obtained  on  U(Y,a)  at  x/OO.l, 


0.2,  0.3  and  0.67  for  6  values  of  a=  0  and  12deg  on  the  top  of  the  figure,  6  and  9deg  during  downstroke  on  the 
left,  and  during  upstroke  on  the  right.  It  can  be  seen  that  the  boundary  layer  velocities  obtained  at  x/C=0.67  are 
not  influenced  all  along  the  period  by  a  bubble  eventually  located  at  this  chord  station.  A  more  accurate  analysis 
of  the  results  shows  that  these  profiles  are  typical  of  a  turbulent  boundary  layer  (U(Y)=kU  <1/n)).  On  the  other 
hand,  velocity  profiles  close  to  the  leading  edge  (x/C=0.1;  0.2;  0.3)  are  affected  by  the  bubble  particularly  at 
x/C=0. 1  when  a=12  and  9deg  during  the  downstroke  and  show  a  lack  of  values  in  the  boundary  layer  (Y<2mm). 
At  such  incidences,  the  bubble  induces  transition  of  the  BL  that  is  characterized  by  scattered  results  inside  the 
B.L.  During  upstroke  the  bubble  effect  weakens  at  x/C=0.1  as  it  can  be  seen  at  a=9deg,  due  to  the  weakness  of 
the  bubble  in  upstroke  at  shown  in  Fig.  13.  Moreover  the  variations  of  the  potential  velocity  (Y>10mm)  with  x/C 
and  incidence  are  in  good  agreement  with  the  variations  of  the  known  pressure  distribution  over  the  upper  side  of 
a  NACA0012  airfoil  (suction  pick  close  to  x/C=0.1).> 

Although  very  few  particles  can  enter  into  the  bubble  (see  Fig.  12),  it  has  been  possible,  allowing  to  a 
long  time  consuming  process,  to  record  velocity  profiles  in  the  bubble.  This  test  was  made  at  x/C=0.2,  a=9deg  in 
downstroke,  where  the  thickness  of  the  bubble  is  2mm  (Fig.  13).  It  is  shown  in  Fig.  15  that  the  velocity  profile 
measured  from  the  wall  to  the  potential  flow  exhibits  a  reverse  flow  due  to  the  vortex  inside  the  bubble.  On  the 
same  figure  is  also  presented  the  velocity  profile  of  a  typical  attached  B.L.,  measured  at  xJ 0=0.3,  where  no 
reverse  flow  is  observed. 


Velocity  profiles  at  high  incidence  (9deg<a<21deg)  and  low  Reynolds  number  (Re=105) 

In  the  case  where  dynamic  stall  occurs  during  the  oscillation  (9deg<a<21deg),  the  separation  region  is 
characterized  by  velocity  records  analogous  to  those  presented  as  long  as  the  bubble  was  concerned  at  low 
incidence.  The  analysis  of  data  lack  on  velocity  realized  at  x/OO.2;  0.3;  and  0.67  allows  determining  at  a  given 
incidence  of  the  oscillation  the  mean  value  of  the  bubble  dimension  normal  to  the  surface  airfoil.  The  dimensions 
of  the  bubble  deduced  from  velocity  records  are  presented  by  white  segments  in  Figure  16  at  incidences  a=12; 
14;  and  18deg.  These  results  compared  to  visualizations  show  a  good  agreement  although  they  concern  averaged 
(velocity  records)  and  instantaneous  (visualizations)  values.  Bubbles  as  high  as  100mm  (Y/C— 0.3)  have  been 
measured. 

Results  obtained  on  the  U(Y,a)  component  at  x/C=0.2;  0.3;  and  0.67  are  presented  in  Figure  17  for  6 
values  of  a: 

-9deg  and  21  deg  on  the  top  of  the  figure,  corresponding  to  maximum  and  minimum  values  of  the 
incidence  during  the  oscillation, 

-15deg  and  18deg  as  typical  examples  of  incidences  during  downstroke  (on  the  left),  and  upstroke  (on 

the  right). 

It  is  worthy  to  note  that  during  downstroke  at  high  incidence  (15deg<a<18deg)  the  separation  can  be 
pointed  out  by  a  no  data  recorded  zone  at  leading  edge  resulting  in  no  value  of  U(Y,a)  on  the  plots  relative  to 
x/OO.2,  and  0.3.  Nevertheless,  a  reverse  flow  due  to  vortices  is  shown  at  x/OO.67,  characterized  by  negative 
values  of  U(Y,a). 

At  minimum  incidence  (9deg),  and  during  upstoke  a=15deg,  it  has  been  possible  to  perform  velocity 
measurements  very  close  to  the  wall;  results  obtained  are  typical  of  a  reattached  flow  with  boundary  layer. 
Profiles  obtained  at  higher  upstroke  incidence  (a=18deg),  reveal  a  bubble  located  at  leading  edge  and  a  flow 
reattachment  at  x/OO.67,  due  to  the  vortex  induced  at  leading  edge  during  upstroke. 


Velocity  profiles  at  high  incidence  (15deg<a<27deg)  and  high  Reynolds  number  (Re— 106) 

Velocity  records  from  laser  velocimeter  performed  at  Re~106  in  S1L  concern  2  chord  stations: 
x/OO.4  and  0.6.  For  each  value  of  x/C,  velocity  measurements  have  been  done  at  the  following 
altitudes: Ymm=80;  60;  40;  20;  10;  8;  6;  4;  2;  1.5;  1;  0.8;  and  0.6. 

Analyzing  the  records  as  those  presented  in  Fig  12  can  carry  out  qualitative  information.  A  region 
where  the  flow  is  reattached  is  characterized  by  a  large  number  of  no  scattered  data  acquired  during  several 
periods  of  oscillation.  The  values  measured  at  the  different  altitudes  allow  to  plot  the  velocity  profiles  of  U,  and 
V,  from  outside  to  the  airfoil  surface.  When  the  wing  is  undergoing  dynamic  stall  with  strong  vorticity  flow 
occurring  on  the  upperside,  records  show  that  the  regions  explored  are  scattered  over  with  less  data.  In  this  case, 
velocity  profiles  become  less  accurate.  Figure  18  presents  as  an  example  data  records  obtained  at  X/C=0.4  (top 
of  the  figure),  and  X/C=0.6  (bottom),  for  Ymm=80.  It  can  be  seen  that  the  region  close  to  21deg  upstroke  is 
scattered  over  with  data,  in  accordance  to  Fig  1 1  where  the  flow  shows  a  high  vorticity  at  this  incidence.  At 
21  deg  downstoke,  the  data  are  gathered  around  the  mean  value  of  the  flow,  attesting  of  a  reattached  flow  as 
shown  in  the  Fig.  11.  Same  comments  can  be  made  at  18deg  up  and  down,  and  16deg  down  (Figs.  8- 10). 
Concerning  27deg,  the  flow  at  Y=80mm  seems  to  be  potential,  as  shown  in  Fig.  1 1  where  Y=80mm  is  very  close 
to  the  boundary  of  the  undisturbed  flow. 


Figure  19  shows  some  examples  of  the  U  (Y,a)  profiles  deduced  from  data  records  as  those  presented 
in  Figl  8  and  concern  incidence  where  the  scatter  of  data  remains  low  enough  between  27deg  and  18deg  in 
downstroke  to  allow  an  accurate  measurement  of  the  velocity.  Fig.  19  gives  the  velocity  profiles  obtained  by 
averaging  the  velocity  results  at  27;  21  and  18deg.  As  the  incidence  decreases  from  21deg  to  18deg  it  becomes 
more  difficult  to  obtain  useful  values  close  to  the  wall,  as  shown  by  the  lack  of  data  on  the  plots. 

CONCLUSIONS 

Two  wings  of  NACA  0012  airfoil  oscillating  in  pitch  has  been  tested  in  2D  and  3D  configurations  at 
Reynolds  Numbers  ranging  from  105  to  106.  The  mean  angles  of  attack  were  a0=6;  15  and  21deg,  the  angular 
amplitude  Aa=±6deg  and  the  frequency  oscillation  k=0.188. 

Smoke  filaments  emitted  from  upstream  visualized  the  flow  on  the  upper  side  of  the  airfoil,  illuminated  by 
a  laser  sheet.  Pictures  of  the  flowfield  were  taken  at  different  instants  of  the  oscillation  by  use  of  a  numerical 
camera  and  movies  were  realized  from  the  pictures. 

Measurements  of  axial  and  radial  components  of  the  velocity  were  performed  using  an  embedded  laser 
velocimeter.  The  velocity  survey  was  made  normally  to  the  airfoil  surface,  from  100mm  outside  to  0.4mm  to 
the  wall. 

The  different  features  of  the  flow  around  the  oscillating  airfoil  has  been  analyzed  both  from  visualizations 
and  velocity  measurements.  It  has  thus  been  possible  to  quantitatively  characterize  a  separation  bubble  close  to 
the  leading  edge  when  the  wing  oscillates  at  low  incidence  (less  than  12deg).  The  bubble  remains  located  at  the 
third  front  part  of  the  upperside  with  a  normal  dimension  to  the  wall  Y/C  varying  between  about  0.02  and  0.01 
during  respectively  downstoke  and  upstroke. 

At  higher  incidence,  pictures  and  velocity  measurements  have  shown  the  different  features  experienced  by 
the  flow  during  downstroke  and  upstroke  as  reattachment,  large  separation  bubble,  and  high  vorticity  flow. 

The  comparisons  made  on  the  results  obtained  in  different  conditions  seem  to  point  out  that  the  effects  of 
Reynolds  numbers  in  the  range  105-106,  and  of  the  tip  wing  are  very  weak,  as  already  remarked  by  W. 
Bouseman*7*. 

Although  velocity  results  presented  in  the  frame  of  the  present  report  concern  only  the  U  component,  a 
database  made  of  both  U  (a,  Y)  and  V  (a,  Y)  has  been  constituted  that  will  be  very  useful  in  computational 
unsteady  boundary  layer. 

This  study  has  clearly  shown  that  embedded  laser  velocimetry  is  a  suitable  and  useful  tool  to  investigate 
unsteady  boundary  layer  and  dynamic  stall  regime  occurring  on  oscillating  airfoils.  In  the  future,  it  is  intended 
to  emphasize  on  the  measurement  of  the  turbulence  components. 
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Steady  case  a=10  deg 


Pitching  case,  incidence  downstroke,  a=10  deg 


Pitching  case,  incidence  upstroke,  a=10  deg 
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Flow  visualization  around  aNACA0012  airfoil,  Re=10^ 


Steady  case  a=\4  deg 


Pitching  case,  incidence  downstroke,  a=14  deg 


Pitching  case,  incidence  upstroke,  a=14  deg 


Flow  visualization  around  aNACA0012  airfoil,  Re=10^ 
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Fig.  6 


Steady  case  oc=10  deg 


Pitching  case,  incidence  downstroke  ,  oc=10  deg 


Pitching  case,  incidence  upstroke  ,  a=10  deg 


Flow  visualization  around  aNACA0012  airfoil,  Re=2  10^ 
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Steady  case  a=14  deg 


Pitching  case,  incidence  downstroke,  a=14  deg 


Pitching  case,  incidence  upstroke,  a=14  deg 
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Flow  visualization  around  aNACA0012  airfoil,  Re=2 
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Steady  case  oc=18  deg 


Pitching  case,  incidence  downstroke,  a=18 
deg 


Pitching  case,  incidence  upstroke  ,  a=18deg 
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visualization  around  aNACA0012  airfoil,  Re= 


Steady  case  a=21  deg 


Pitching  case,  incidence  downstroke,  a=21  deg 


Pitching  case,  incidence  upstroke,  a=21deg 
Flow  visualization  around  aNACA0012  airfoil,  Re 


Steady  case  a=16  deg 


Pitching  case,  incidence  downstoke,  a=16deg 


Flow  visualization  around  aNACA0012  airfoil,  Re-  10 
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Steady  case  oc=27  deg 


Pitching  case,  a=27  deg 
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Flow  visualization  around  aNACA0012  airfoil,  Re—  10^ 
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Bubble  dimension  with  incidence 

ctQ=6;  Aa=±6;  Re=10^;  k=0.188 
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SEPARATION  BUBBLE:  VISUALIZATION  AND 
THICKNESS  DEDUCED  FROM  VELOCITY  RECORDS 


a0=  16deg;  Aa=6  deg;  k=0.188;  Re=105 
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